Coastal  Engineering  57  (2010)  610-619 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Coastal  Engineering 

journal  homepage:  www.elsevier.com/locate/coastaleng 


a 

Coastal 

Engineering 


Calculation  of  beach  change  under  interacting  cross-shore  and  longshore  processes 

Hans  Hanson a  *.  Magnus  Larson a,  Nicholas  C.  Kraus  b 

a  Water  Resources  Engineering,  Lund  University,  Box  118,  22100  Lund,  Sweden 

b  U.S.  Army  Engineer  Research  and  Development  Center,  Coastal  and  Hydraulics  Laboratory,  3909  Halls  Ferry  Road,  Vicksburg,  MS  39180-6199,  United  States 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  13  July  2009 

Received  in  revised  form  18  January  2010 

Accepted  4  February  2010 

Available  online  7  March  2010 


Keywords: 

Longshore  sediment  transport 

Beach  response 

Groins 

Coastal  structures 
Numerical  modeling 
Shoreline  evolution 
Dune  erosion 
Windblown  sand 


This  paper  presents  a  mathematical  approach  and  numerical  model  that  simulates  beach  and  dune  change  in 
response  to  cross-shore  processes  of  dune  growth  by  wind  and  dune  erosion  by  storms,  and  by  gradients  in 
longshore  sand  transport  that  will  alter  shoreline  position.  Sub-aerial  transport  processes  are  represented, 
whereas  sub-aqueous  transport  is  neglected.  The  system  is  tightly  coupled  morphologically,  with  the  berm 
playing  a  central  role.  For  example,  the  potential  for  sand  to  be  transported  to  the  dune  by  wind  depends  on 
berm  width,  and  sand  lost  in  erosion  of  the  dune  during  storms  can  widen  the  berm.  Morphologic 
equilibrium  considerations  are  introduced  to  improve  reliability  of  predictions  and  stability  of  the  non-linear 
model.  An  analytical  solution  is  given  under  simplification  to  illustrate  properties  of  the  model.  Sensitivity 
tests  with  the  numerical  solution  of  the  coupled  equations  demonstrate  model  performance,  with  one  test 
exploring  beach  and  dune  response  to  potential  increase  in  storm-wave  height  with  global  warming.  Finally, 
the  numerical  model  is  applied  to  examine  the  consequences  of  groin  shortening  at  Westhampton  Beach, 
Long  Island,  New  York,  as  an  alternative  for  providing  a  sand  supply  to  the  down-drift  beach.  Results  indicate 
that  the  sand  will  be  released  over  several  decades  as  the  shoreline  and  dune  move  landward  in  adjustment 
to  the  new  equilibrium  condition  with  the  shortened  groins. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Coastal  management  requires  quantitative  predictive  capability 
that  represents  episodic  erosion  by  storms,  dune  build-up  by  wind, 
and  long-term  change  produced  by  gradients  in  longshore  sand 
transport.  Regional  processes  such  as  shadowing  by  adjacent  large 
land  masses  and  sand  storage  and  transfer  at  inlets  also  exert  control 
of  beach  change  at  the  project  scale  (Larson  et  al,  2002a, b).  Realistic 
representation  of  such  processes  therefore  requires  their  mathemat¬ 
ical  descriptions  and  associated  numerical  modeling  on  several  spatial 
and  temporal  scales.  The  interaction  between  waves,  structures,  and 
morphological  processes  must  be  described,  but  within  a  scheme  that 
is  compatible  with  boundary  conditions  and  processes  at  regional 
scale.  Also,  model  robustness  and  reliability  must  be  assured  by 
introducing  equilibrium  concepts  towards  which  evolution  is  striving 
rather  than  using  an  open-ended  modeling  system  of  non-linear 
coupled  processes. 

This  paper  presents  a  theory  for  calculating  beach  and  dune 
change  through  coupled  storm  and  longer  term  or  typical  longshore 
processes.  Change  in  the  location  of  the  dune  toe,  berm  width,  and 
shoreline  position  are  calculated,  while  maintaining  longshore 
transport  rates  representative  of  the  regional  long-term  behavior  of 
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the  coast.  The  paper  begins  with  an  introduction  to  and  representa¬ 
tion  of  dune  evolution,  which  is  then  connected  to  longshore 
transport  processes  as  represented  in  one-line  models  such  as 
GENESIS  (Hanson,  1989;  Hanson  and  Kraus,  1989). 

2.  Representation  of  dune  evolution  processes 

Dunes  are  formed  by  wind-blown  sand  transport  from  the  berm 
and  foreshore.  Psuty  (1990)  discusses  dune  transformation  processes 
for  the  general  area  of  Long  Island,  NY,  that  enters  the  case  study 
below.  The  berm  refers  to  the  mildly  seaward  sloping,  nearly- 
horizontal  portion  of  the  sub-aerial  beach  seaward  of  the  dune.  The 
portion  of  the  beach  regularly  exposed  to  waves  in  the  swash  zone, 
sloping  downwards  towards  the  still-water  shoreline,  is  referred  to  as 
the  foreshore.  Fig.  1  is  a  definition  sketch  introducing  beach 
morphology  terminology  and  notation  employed.  In  the  present 
work,  the  subaqueous  part  of  the  profile  is  described  by  an 
equilibrium  shape  according  to  Dean  (1977).  Gradients  in  longshore 
transport  are  assumed  to  act  upon  the  profile  from  berm  crest  to 
depth  of  closure  (location  denoted  byyc). 

Waves  approaching  the  beach  during  high  water  level  may  attack 
the  dune  during  storms,  causing  erosion  and  recession  of  the  dune 
line.  Sand  released  from  the  eroded  dune  supplies  the  beach  with 
material  that  can  be  deposited  offshore,  especially  during  long- 
duration  storms.  Larson  and  Kraus  (1991)  examined  storm-induced 
beach  change  for  differing  durations  of  storm  surge  with  the  SBEACH 
model  (Larson  and  Kraus,  1989)  and  found  that  dune  erosion  by  wave 
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attack  depend  on  berm  width,  with  wider  berm  provides  increased 
protection  for  the  dune  against  wave  attack.  Build-up  of  dunes  by 
wind-blown  sand  transport  also  depends  on  berm  width  (Hotta, 
1984),  for  which  a  larger  width  promotes  transport  by  wind.  Hotta 
(1984)  showed  that,  in  his  studies,  a  certain  distance,  about  10  m,  was 
required  for  transport  by  wind  to  become  fully  developed.  This  section 
describes  a  new  theory  that  allows  for  numerical  modeling  of  the 
interaction  between  the  dune  and  the  berm,  including  consequences 
for  shoreline  evolution. 

2.1.  Dune  build-up  by  wind 

It  is  assumed  that  a  dune  grows  by  sand  blown  landward  from  the 
berm,  which,  in  turn,  exchanges  material  with  the  foreshore  and  the 
underwater  profile.  A  one-dimensional  sand  conservation  equation 
for  the  dune  volume  growth  by  windblown  transport,  assuming  that 
all  the  blown  sand  is  trapped  by  the  dune,  is: 


dVow 

dt 


Qw 


(1) 


beaches,  a  stronger  increase  for  wider  beaches,  and  an  upper  limit  that 
is  approached  gradually  for  wide  beaches,  while  at  the  same  time 
providing  a  continuous  description  of  the  transport  with  changes  in 
berm  width,  is, 


1-0.5  1-tanh 


Qgrad 


(Tb  Vd  y50) 


(4) 


where  qwo  is  the  maximum  transport  by  wind  for  an  infinitely  wide 
beach,  dependent  on  water  and  sand  properties,  yB  and  yD  are  the 
distances  to  the  seaward  end  of  the  berm  and  the  dune  toe, 
respectively  (see  Fig.  1),  with  the  y-axis  pointing  offshore,  y50  is  the 
distance  from  the  seaward  end  of  the  berm  to  where  the  wind-blown 
transport  has  reached  50%  of  its  maximum,  and  qgmd  is  the  transport 
gradient  at  y50.  Bagnold  (1954)  suggested  the  transport  rate 
relationship  qwo  =  I<wu*/g,  where  u*  is  the  wind  shear  velocity,  g  is 
the  acceleration  due  to  gravity,  and  I<w  is  an  empirical  coefficient  that 
quantifies  the  influence  of  sand  properties  on  the  transport  rate. 
Eqs.  (4)  and  (3)  describe  a  dune  that  advances  towards  the  berm  crest, 
although  the  rate  of  advance  will  decrease  with  time  as  the  berm 
width  decreases 


where  VDw  is  the  dune  volume  (subscript  w  refers  wind),  t  is  time,  and 
qw  is  the  onshore  sand  transport  by  wind.  If  the  dune  height  Dd  above 
the  berm  elevation  remains  constant  during  growth  together  with  the 
dune  shape,  the  seaward  movement  of  the  dune  toe  AyDw  for  a  given 
increase  in  dune  volume  AVDw  is: 


2.2.  Erosion  by  wave  impact 

During  storms,  waves  may  strike  and  erode  the  dune  during  high 
water  level.  Denoting  the  erosion  of  the  dune  due  to  wave  impact  as  q0,  the 
associated  dune  volume  change  (loss)  AVDo  during  a  small  time  At  yields: 


AyDw  = 


AVdw 

Dd 


(2) 


(5) 


Expressing  Eq.  (1)  in  terms  of  dune  toe  advance  yields: 

dyDw  _  Qw  n\ 

dt  -  Dd 

It  is  assumed  that  sand  transport  to  the  dune  is  related  to  the  width 
of  the  berm  up  to  some  distance  over  which  equilibrium  conditions 
have  developed,  implying  that  beyond  equilibrium  a  wider  beach  does 
not  generate  more  transport  by  wind,  (Davidson-Arnott  and  Law, 
1990,  Davidson-Arnott  et  al.,  2005;  see  recent  review  of  wind-blown 
sand  transport  prediction  by  Hopf  and  Sherman,  2007).  A  simple 
equation  that  exhibits  these  properties,  with  a  slow  but  gradual 
increase  of  the  transport  rate  by  wind  with  beach  width  for  narrow 


It  is  assumed  that  no  overwash  occurs  and  that  the  dune  is  not 
completely  eroded  (i.e.,  no  breaching).  The  erosion  rate  due  to  dune 
impact  by  waves  may  be  estimated  as  (Larson  et  al.,  2004): 


Qo  =  4Q 


(R  +  A  h-zD)2 

T 


R>zD-Ah 


(6) 


where  R  is  the  run-up  height  (including  setup)  estimated  from 
R  =  aVH0L0l  in  which  H0  is  the  deepwater  root-mean-square  wave 
height,  L0  is  the  deepwater  wavelength,  and  a  is  a  coefficient  (about 
0.15,  which  corresponds  to  a  representative  foreshore  slope);  Ah  is 
the  surge  level  (including  tide  elevation  relative  to  mean  sea  level 
(MSL));  zD  is  the  dune  toe  elevation  (with  respect  to  MSL);  T  is  the 
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swash  period  (taken  to  be  the  same  as  the  wave  period);  and  Cs  is  an 
empirical  coefficient.  In  the  numerical  implementation,  q0  varies  at 
each  time  step  and  is  computed  following  Eq.  (6)  from  the  input  time 
series  of  waves. 

Eq.  (6)  is  derived  based  on  the  wave  impact  on  the  dune,  where 
the  speed  of  the  wave  striking  the  dune  is  a  key  parameter  (estimated 
through  R).  If  there  is  a  wide  berm  in  front  of  the  dune,  the  up-rushing 
wave  may  lose  energy  as  it  travels  across  the  berm  and  strike  the  dune 
with  less  speed,  causing  less  impact  and  erosion.  Such  an  energy  loss 
is  not  included  in  Eq.  (6).  In  this  study,  a  ballistics  approach  was  taken 
to  describe  the  front  speed  of  the  wave  affected  by  friction  as  it 
propagates  towards  the  dune  face  (Erikson  et  al.,  2005).  The 
governing  equation  for  the  wave  front  speed  is: 


transport  by  wind  (constructive)  and  wave  impact  (destructive  to  the 
dune): 


dVD  dVDw  dVDo 
dt  dt  dt 


fdypw  ,  typo) 
(  dt  dt  J 


Do  =  ^ Dd  =  qB 


Thus,  if  qw>q0  the  dune  will  advance,  whereas  if  qw<q0  it  will 
recede,  and  at  qw  =  q0  equilibrium  occurs. 

If  a  volume  AVD  is  added  to  the  dune,  the  same  volume  should  be 
taken  from  the  profile  over  its  active  depth,  which  is  the  sum  of  the 
vertical  distance  from  the  berm  crest  DB  to  the  still-water  level  and  the 
depth  of  closure  Dc  (see  Fig.  1 ).  Under  the  assumption  that  dune  and 
beach  profile  change  occur  with  the  same  shape,  continuity  requires 
that: 


AyDUD  +  AyB{DB  +  Dc)  —  0 


(12) 


where  u0  is  the  front  speed  and  kf  a  friction  coefficient.  Rewriting 
Eq.  (7)  in  terms  of  a  coordinate  s  along  the  berm,  originating  at  the 
berm  crest  and  pointing  onshore,  yields  (note  that  du0/dt=  u0dii0/ds): 


where  A yB  is  the  berm  crest  translation  due  to  interaction  with  the 
dune  system.  Eq.  (12)  may  be  re-arranged  to  yield: 


u 


0 


duo 

ds 


-kfu 


iuo 


(8) 


With  the  boundary  condition  u0  =  uc  for  s  =  0,  the  solution  to 
Eq.  (8)  is, 


u0  =  uc  exp  (—kfs'j  (9) 

where  uc  is  the  velocity  of  the  up-rushing  wave  at  the  berm  crest 
(which  is  equivalent  to  the  velocity  at  the  dune  toe,  if  friction  is 
neglected).  The  effect  of  this  velocity  reduction  on  q0  is  quantified  by 
an  adjusted  run-up  height  R'  in  Eq.  (6)  instead  of  R,  where  R ’  is 
calculated  from: 


R'  =  R  exp  2fysB)  +  zD[l  — exp(— 2fysB)]  (10) 

where  sB=yB—yD.  This  expression  gives  a  more  realistic  description 
than  Eq.  (6),  because  it  accounts  for  the  reduction  in  impact  as  a  result 
of  horizontal  travel  distance  of  the  wave  front. 

If  /</— ►  0,  R'  — >  R,  and  there  is  no  effect  of  friction  on  the  run-up 
height  and  the  impact  velocity.  On  the  other  hand,  if  R'-»zD 

and  the  run-up  only  reaches  the  berm  crest,  giving  an  impact  velocity 
of  zero.  The  friction  coefficient  will  be  a  function  of  the  flow  properties 
and  the  roughness  of  the  sand  surface  on  the  berm  (Camenen  et  al., 
2006,  2009),  and  will  in  the  general  case  vary  across  the  berm. 
However,  for  the  schematic  approach  taken  here,  the  friction  coef¬ 
ficient  is  specified  as  a  constant  value.  Studies  on  bore  propagation 
on  low-sloping  foreshores  may  give  insight  to  improved  values  on 
the  friction  coefficient.  Raubenheimer  et  al.  (1995)  used  a  coefficient 
value  of  /=  0.008  in  solving  the  shallow-water  equations  in  the 
swash.  They  stated  that  model  simulations  were  not  sensitive  to 
values  in  the  range  of  0.005  </<  0.025.  It  is  difficult  to  compare 
/  with  kf,  because  kf  is  dimensional  in  Eq.  (7)  and  needs  to  be 
divided  with  the  local  water  depth  to  be  consistent  with  /.  In  the 
present  study  a  value  of  kf=  0.02  was  employed,  which  for  a  realistic 
bore  height  puts  the  value  of  /  in  the  lower  range  of  previous 
investigations. 

2.3.  Coupling  between  subaerial  profile  response  and  shoreline  change 

The  resulting  volume  change  of  the  dune  AVD  and  the  associated 
movement  of  the  location  of  the  dune  toe  A yD  are  the  combination  of 


AyB  =  -AyDp—p 


(13) 


This  equation  provides  a  simple  estimate  of  the  needed  profile 
recession  due  to  cross-shore  processes,  from  the  foot  of  the  dune  to 
the  depth  of  closure,  to  produce  a  certain  dune  advance,  and  vice 
versa. 

Next,  the  cross-shore  exchange  between  the  berm  and  dune  is 
combined  with  the  alongshore  sand  transport  rate  caused  by 
obliquely  breaking  waves  through  the  continuity  equation  of 
shoreline  change: 


dy  =  1  fdQ)  dyg 

dt  Db  +  Dc  (dx)  dt 


1 

Db  +  Dc 


dx  q° 


(14) 


where  the  berm  translation  due  to  cross-shore  interaction  between 
the  dune  and  berm  are  linearly  added  to  the  contribution  by  the 
gradient  in  longshore  transport  rate,  dQJdx ,  to  obtain  the  total 
shoreline  response. 


3.  Analytical  solution  to  equations  governing  dune  processes 


Based  on  the  theoretical  developments  and  empirical  relationships 
discussed  in  the  previous  section,  the  governing  equations  for  the 
movement  of  the  dune  toe  yD  and  the  berm  crest  yB  location  are: 

_  Qw  Qo  /-|  r\ 

dt  -  Dd 


dy b  Qw  Qo 

dt  Db  +  Dc 


(16) 


If  the  expressions  for  dune  growth  by  wind-blown  sand  transport 
and  by  erosion  caused  by  wave  impact  are  simple,  analytical  solutions 
may  be  developed  to  solve  Eqs.  (15)  and  (16).  Such  an  approach 
provides  insight  to  the  general  properties  of  the  governing  equations 
and  the  time  scale  controlling  the  response  of  the  dune  and  its 
coupling  to  the  beach.  It  is  not  possible  to  find  analytical  solutions  if 
the  wind-blow  sand  transport  rate  is  described  by  Eq.  (4).  Thus,  a 
simpler  form  is  introduced  to  proceed  with  an  analytical  solution.  An 
equation  that  exhibits  similar  properties  to  Eq.  (4),  that  is,  growth  in 
the  transport  rate  with  distance  from  the  berm  crest  towards  a 
maximum  value,  is: 

Qw  =  <Jwo[1-exp(-a(yB-yD))] 


(17) 
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where  a  is  a  scale  coefficient.  Erosion  by  wave  impact  is  described  by 
Eq.  (6),  assuming  that  the  wave  and  water  level  conditions  are 
constant  and  that  energy  losses  are  negligible  along  the  berm. 

Combining  Eqs.  (15)  and  (16),  substituting  in  Eq.  (17)  for  qw,  and 
introducing  the  variable  w=yB— yD  (representing  the  berm  width) 
produces  the  following  equation  to  solve: 


dw 

~dt 


^(1  +  B)[A- 

UD 


■  exp  (—aw)] 


where  A  and  B  are  defined  by: 


_  Qwo  Qo  _  |  _  do 


Do  +  Dr 


(18) 

(19) 

(20) 


Fig.  2.  Non-dimensional  evolution  of  the  berm  width  for  varying  q0/qwo-ratios  in  the 
In  deriving  the  solution  to  Eq.  (18),  it  is  convenient  to  define  the  case  of  dune  build-up  by  wind  and  dune  erosion  due  to  wave  impact, 
following  non-dimensional  variable: 


£,  =  exp[-a(yB-yD)]  =  exp  (-aw) 

The  solution  to  Eq.  (18)  may  be  written: 

£  =  _ _ 

(A— yexp(— cM(l  +  B)qwot  /  Dd)  + 

where: 


£,„  =  exp  [— a(yBo— yDo)]  =  exp(-aw0) 


(21) 


(22) 


(23) 


and  subscript  o  denotes  conditions  at  t  =  0. 

If  t—>° o  in  Eq.  (22),  > A  for  A>0,  whereas  for  A<0  there  is  no 

asymptotic  solution,  and  erosion  continues  indefinitely.  A  positive 
value  on  A  implies  that  the  potential  transport  rate  by  wind  is  larger 
than  the  erosion  rate  due  to  wave  impact  (qW0>q0‘,  see  Eq.  (19)), 
yielding  an  equilibrium  width  we  of  the  berm  given  by: 


=  (yB-yD). 


In 


f  __Qwo_\ 
\Qwo  QoJ 


(24) 


where  the  subscript  e  denotes  equilibrium.  This  expression  can  also  be 
obtained  by  setting  dw/dt  =  0  in  Eq.  (18). 

To  display  the  properties  of  Eq.  (22),  the  equation  is  developed  and 
written  in  non-dimensional  form  according  to: 


w'  =  ^ln  Q[l— (1—  A  exp  (a'))  exp  (— a'ABt')]^ 


(25) 


where  w'  =  w/w0  =  (yB-yD)/(yBo-yDo),  t' =  qwot/w0DD,  and  a' = 
aw0.  Fig.  2  illustrates  the  time  evolution  of  the  berm  width  towards 
equilibrium  for  various  values  on  the  ratio  q0/qwo  ami  a' =  1.0. 
The  value  8  =  1/3  represents  a  dune  height  Dd  that  is  one-third  that 
of  (Db  +  Dc).  The  initial  berm  width  is  shorter  than  the  equilibrium 
value  for  the  largest  q0/qwo-ratio,  which  causes  the  berm  width 
to  increase.  Fig.  3  displays  calculated  evolution  of  the  berm,  where 
a' =  1.0  and  A  =  0.6  (corresponding  to  qo/qwo  =  0A).  A  relatively 
larger  active  profile  height  with  respect  to  the  dune  height  causes  a 
slower  response  towards  equilibrium. 

With  analytical  solutions  to  somewhat  simplified  versions  of  the 
coupled  equations  demonstrating  reasonable  results,  the  governing 
equations  for  coupled  cross-shore  and  longshore  processes  are  now 
solved  numerically  for  a  field  application. 


4.  Application  to  Westhampton  Beach,  Long  Island,  New  York 

Westhampton  is  located  on  the  eastern  south  shore  of  Long  Island, 
between  Shinnecock  Inlet  and  Moriches  Inlet  (Fig.  4).  The  beach  is 
composed  of  medium  sand  with  a  representative  median  grain  size 
d50  of  0.4  mm,  and  the  net  average  annual  rates  of  longshore  sand 
transport  at  the  site  have  been  variously  estimated  at  115,000- 
230,000  m3/yr  directed  to  the  west,  although  reversals  are  typical  in 
the  summer  and  may  persist  through  a  particular  year  (Rosati  et  al., 
1999,  Larson  et  al.,  2002a).  Long-term  wave  and  water  level  data  were 
input  to  the  numerical  model  to  calculate  dune  change,  berm  change, 
and  shorelines  position.  Hindcast  data  from  the  U.S.  Army  Corps  of 
Engineers  Wave  Information  Study  (WIS)  (http://chl.erdc.usace.army. 
mil/wis)  were  available  from  01/01/76  to  31/12/95  and  measured 
water  levels  from  01/01/80  to  31/12/00,  both  at  3-hour  intervals. 
Thus,  the  period  with  simultaneous  wave  and  water  level  data 
covered  the  period  01/01/80  to  31/12/95.  Based  on  examination  of  a 
large  number  of  profile  surveys  along  the  beach,  the  average  dune 
height  (above  the  berm)  was  estimated  to  be  1.5  m,  the  average  berm 
elevation  3  m  above  MSL,  and  the  average  depth  of  closure  8  m.  The 
transport  rate  by  wind  at  equilibrium,  qwo,  was  set  to  0.0025  kg/m/s, 
based  on  analysis  of  the  actual  dune  toe  migration  on  site.  This 
magnitude  is  in  good  agreement  with  the  rates  presented  in 
Davidson-Arnott  and  Law  (1990),  who  measured  transport  by  wind 


Fig.  3.  Non-dimensional  evolution  of  the  berm  width  for  varying  6-values  in  the  case  of 
dune  build-up  by  wind  and  dune  erosion  due  to  wave  impact. 
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Fig.  4.  Long  Island,  NY,  location  map  and  estimated  net  longshore  sand  transport  rates. 


across  the  berm  in  the  field.  This  transport  rate  corresponds  to  a  wind 
speed  of  slightly  more  than  6  m/s. 

In  the  following,  two  examples  based  on  the  general  conditions  at 
Westhampton  are  first  given,  neglecting  the  groin  field  there,  to 
illustrate  behavior  and  trends  in  the  numerical  solution.  An 
application  of  the  model  is  then  presented  where  groin  shortening 
is  considered  as  an  engineering  alternative  for  supplying  sand  to  the 
down-drift  (western)  beach. 


Eq.  (13).  The  calculated  shoreline  variation  is  about  6  m  over  the 
simulation  period,  which  is  significantly  smaller  than  values  observed, 
probably  because  wave-induced  erosion  of  the  berm  is  not  considered 
represented. 

In  a  second  configuration,  the  shoreline  (thick  solid  line)  was  set 
back  30  m,  compared  to  the  initial  conditions  in  the  first  configura¬ 
tion,  while  leaving  the  initial  dune  toe  location  the  same  as  in  the  first 
configuration.  The  dune  (thick  dashed  line)  displays  similar  short¬ 
term  fluctuations  as  in  the  first  configuration,  but  over  the  years  it 
gradually  recedes  slightly  more  and  recovers  slightly  less  compared  to 
the  first  configuration,  attaining  its  equilibrium  position  further  back. 
The  difference  between  the  dune  toe  locations  for  the  two  configura¬ 
tions  is  indicated  by  the  thick,  dash-dotted  line  in  Fig.  5.  Initially,  the 
difference  is  zero  but  gradually  it  increases  over  the  years.  After  some 
22  years,  the  difference  has  grown  to  26  m,  after  which  it  remains 
virtually  constant.  Over  the  same  time,  the  difference  in  shoreline 
position  (thin,  dash-dotted  line)  decreased  from  its  original  30  m  to 
the  same  26  m  as  for  the  dune.  Thus,  the  beach  and  dune  were  re¬ 
established  in  a  dynamic  equilibrium  where  the  beach  and  dune  have 
the  same  relative  locations  as  before  the  set  back  of  the  shoreline. 
These  tests  demonstrate  that  the  model  of  the  cross-shore  exchange 
between  the  beach  and  the  dune  produces  a  long-term  dynamic 
equilibrium  balance.  If  the  berm  width  is  perturbed  from  its 
equilibrium  value,  the  beach  and  dune  will  change  accordingly  to 
re-establish  a  new  dynamic  equilibrium. 


4.1.  Beach  change  by  cross-shore  processes  only 

To  investigate  model  coupling  between  the  dune  and  the  shoreline 
with  only  cross-shore  processes  active,  several  sensitivity  tests  were 
run  for  a  30-year  period  specifying  wave  uniformity  alongshore.  This 
configuration  allows  examination  of  predicted  time-varying  shoreline 
and  dune  foot  locations  in  an  arbitrary  beach  profile.  The  shoreline 
(thin  solid  line  in  Fig.  5)  was  initially  located  360  m  from  the  baseline 
and  the  dune  toe  (thin  dashed  line)  280  m  from  the  same  baseline, 
resulting  in  an  initial  berm  width  of  80  m.  This  initial  berm  width  was 
close  to  the  value  corresponding  to  its  long-term  equilibrium,  as 
indicated  by  the  absence  of  a  clear  long-term  trend  in  the  shoreline 
position,  although  small  shoreline  fluctuations  occur  at  shorter  time 
scales.  Similarly,  the  location  of  the  dune  toe  shows  no  long-term 
trend.  However,  location  of  the  dune  toe  varies  more  than  the 
shoreline  position  because  the  active  beach  profile  height  (8  +  3  m)  is 
more  than  seven  times  larger  than  the  dune  height  (1.5  m),  c.f. 


4.2.  Beach  change  by  combined  cross-shore  and  longshore  processes 

To  illustrate  the  interaction  between  longshore  and  cross-shore 
transport  processes,  an  example  is  given  for  a  straight  shoreline  and 
longshore  transport  generated  using  the  same  WIS  wave  hindcast 
time  series  for  Westhampton  Beach  as  in  the  previous  example.  The 
simulation  covers  the  15-year  period  July  1, 1980  to  June  30, 1995.  A 
constant  gradient  was  imposed  on  the  longshore  transport  such  that 
the  beach  would  erode  when  the  transport  direction  was  positive  (to 
the  west)  and  accrete  for  negative  transport  rates.  The  simulation 
results  are  illustrated  in  Fig.  6,  where  the  shoreline  evolution  with 
only  longshore  processes  active  is  shown  as  a  solid  line.  The  calculated 
net  shoreline  change  was  7.5  m  corresponding  to  0.5  m/yr,  which  is  a 
realistic  number  for  many  eroding  sites  along  the  eastern  part  of  the 
south  shore  of  Long  Island  (Rosati  et  al.,  1999).  Because  cross-shore 
transport  processes  were  not  included  in  this  simulation,  the  location 
of  the  dune  toe  did  not  move. 
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Fig.  5.  Temporal  evolution  of  shoreline  and  dune  toe  location  for  two  sensitivity  tests  encompassing  cross-shore  processes  only. 
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Fig.  6.  Shoreline  change  over  15  years  with  a)  only  longshore  (LS)  processes  included  (solid  line),  and  b)  both  longshore  and  cross-shore  (CS)  processes  included  (dashed  and  dash- 
dotted  lines). 


If  cross-shore  transport  is  included  through  wind-blow  sand 
transport  and  dune  erosion,  dune  volume  will  vary  which,  in  turn,  will 
induce  increasing  shoreline  fluctuations  as  material  exchange  occurs 
between  the  dune  and  the  shoreline  (Fig.  6,  dashed  line).  The  long¬ 
term  longshore  trend  is  still  apparent  with  a  shoreline  net  change 
trend  of  6.7  m  over  the  15  years.  As  seen  from  Eq.  (13),  shoreline  and 
dune  toe  fluctuations  are  coupled.  For  Westhampton  Beach,  dune 
height  Dd  is  on  the  average  1.5  m,  the  berm  height  DB  is  3  m,  and  the 
depth  of  closure  Dc  is  8  m.  Thus,  dune  fluctuation  will  be  scaled  as  1 1/ 
1.5  times  the  shoreline  fluctuations,  as  seen  from  the  Fig.  6. 

To  investigate  the  role  of  wave  sequencing  and  wave  variability  on 
shoreline  change,  the  same  schematic  case  as  above  was  run  for  a 
series  of  different  wave  climates.  In  addition  to  the  original  WIS 
hindcast  time  series,  ten  new  1 5-year  time  series  were  constructed  by 
random  sampling  of  individual  years  from  the  original  series  (by 
sampling  complete  years;  seasonality,  typical  annual  storm  waves, 
and  associated  storm  water  level  were  maintained).  Sampling  was 
done  with  replacement,  so  that  a  specific  individual  year  could  be 


sampled  more  than  once.  Fig.  7  displays  calculated  shoreline  evolution 
for  only  longshore  processes.  The  thick  line  represents  an  average  of 
the  1 1  individual  shoreline  locations.  For  this  average  shoreline,  the 
net  shoreline  change  over  the  15  years  was  7.3  m,  i.e.,  close  to  the 
change  for  the  original  wave  time  series.  The  difference  in  shoreline 
location  after  1 5  years  is  2.3  m  between  the  run  with  most  net  erosion 
and  the  one  with  the  least  net  erosion,  corresponding  to  32%  of  the 
average  change. 

Fig.  8  shows  the  shoreline  evolution  for  the  same  wave  time  series, 
but  with  cross-shore  processes  included.  There  is  considerably  more 
variability  in  shoreline  location  for  each  wave  climate  as  well  as 
between  wave  climates,  because  of  cross-shore  exchange  of  material 
between  the  dune  and  the  shoreline.  For  the  average  shoreline,  the 
net  shoreline  change  over  the  1 5  years  was  5.9  m,  slightly  less  than  for 
longshore  processes  alone,  because  dune  erosion  supplies  more  sand 
to  the  berm,  than  the  longshore  gradient  in  transport  is  able  to 
remove.  At  the  same  time,  the  difference  in  shoreline  location  after 
1 5  years  is  as  much  as  2.4  m  between  the  run  with  most  net  erosion 


Fig.  7.  Shoreline  change  caused  by  longshore  processes  only  over  15  years  for  1 1  different  wave  climates.  The  thick  line  represents  an  average  of  the  1 1  individual  shoreline  locations. 
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Fig.  8.  Shoreline  change  caused  by  combined  longshore  and  cross-shore  processes  over  15  years  for  11  different  wave  climates.  The  thick  line  represents  an  average  of  the  11 
individual  shoreline  locations. 


and  the  one  with  the  least  net  erosion,  corresponding  to  41%  of  the 
average  change. 

From  a  coastal  management  point  of  view,  variability  in  shoreline 
location  about  some  mean  must  be  considered  (Smith  and  Kraus, 
1999).  For  determining  the  location  of  housing  and  infrastructure 
development,  the  shoreward-most  position  of  the  shoreline  over 
some  period  of  time  with  some  suitable  set  back  might  be  of  primary 
interest,  whereas  in  considering  potential  bypassing  around  coastal 
inlets  and  their  jetties,  the  seaward-most  location  might  be  of  interest. 
Thus,  to  represent  these  changes  on  different  time  scales,  long-term 
and  short-term  processes  must  be  combined  in  a  modeling  system, 
where  the  former  is  more  closely  related  to  longshore  transport  and 
transport  by  wind,  and  the  latter  to  episodic  cross-shore  transport 
that  may  occur  during  storms  and  elevated  water  level. 

Another  issue  with  relevance  for  global  climate  change  is  the 
sensitivity  of  shoreline  evolution  to  changes  in  wave  height.  To 
investigate  how  the  model  responds  to  variations  in  wave  height,  runs 
were  made  over  15  years  with  seven  different  wave  climates,  keeping 


the  water  level  constant.  In  addition  to  the  original  wave  time  series, 
six  new  time  series  were  constructed  by  increasing  the  offshore  wave 
height  10, 20,  and  30%,  respectively,  and  by  reducing  the  wave  height 
with  the  same  percentage  (Fig.  9).  After  about  1  year,  a  storm  with 
larger  waves  occurs,  inducing  dune  erosion  that  produces  significant 
shoreline  advance.  Longshore  transport  gradients  cannot  keep  pace 
with  (remove)  the  supply  of  sand  from  the  dunes  initially,  but 
gradually  longshore  transport  erodes  the  beach.  In  fact,  the  erosion 
caused  by  longshore  transport  will  eventually  be  larger  for  the  higher 
waves.  Thus,  the  effect  of  increased  sand  supply  from  the  dunes  to  the 
berm  system  will  disappears  after  about  9  years. 

Subsequently,  the  beach  is  struck  by  another  large  storm. 
However,  this  time  the  wave  climate  with  the  smaller  waves  (dash- 
dot  lines)  will  induce  greater  change  than  the  more  severe  wave 
climate  (dashed  lines)  indicated  by  a  larger  vertical  displacement  of 
the  dash-dot  lines.  The  explanation  is  that  due  to  larger  sand  losses 
from  the  dune  by  greater  cross-shore  transport  for  the  larger  waves, 
the  system  contains  less  sand  available  for  transport  in  the  large  wave 


Fig.  9.  Shoreline  change  caused  by  combined  longshore  and  cross-shore  processes.  Runs  were  made  over  1 5  years  using  7  different  wave  climates.  In  addition  to  the  original  WIS  time 
series,  6  new  time  series  were  constructed  by  increasing  the  offshore  wave  height  10,  20,  and  30%  (thick  lines),  respectively  and  by  reducing  the  wave  height  with  the  same 
percentage  (thin  line). 
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cases.  Thus,  even  though  the  shoreline  locations  coincide,  the  dune 
toe  location  is  further  landward  for  the  larger  waves  and  higher  water 
levels.  In  spite  of  the  waves  and  the  associated  run-up  being  larger, 
wave  impact  on  the  dune  becomes  smaller  from  the  larger  waves 
because  the  dune  toe  is  located  further  landward.  This  response  of  the 
dune  and  the  shoreline  illustrates  the  complexity  of  the  system  and 
the  significant  role  of  sequencing  regarding  the  forcing. 

4.3.  Evaluation  of  groin  shortening 

There  are  presently  16  groins  at  Westhampton  Beach.  The  first 
eleven  and  easternmost  (up-drift)  groins  were  constructed  over 
1965-1966  in  response  to  breaching  during  the  March  1962 
northeaster.  Four  additional  groins  were  constructed  to  the  west 
during  1969-1970.  The  rubble  stone  groins  in  this  field  of  15 
structures  were  nominally  146  m  long,  tapering  to  the  sea  bottom 
and  having  spacing  of  approximately  365  m.  Kraus  and  Batten  (2006) 
provide  the  most  recent  summary  of  groin  configuration  and 
functioning  at  Westhampton  Beach. 

The  original  groins  functioned  as  intended  (Nersesian  et  al.,  1992) 
in  protecting  a  once-vulnerable  5-km  long  segment  of  barrier  beach 
that  had  experienced  breaching  during  the  September  1938  Great 
New  England  Hurricane,  a  1958  storm,  and  the  March  1962 
northeaster.  Although  the  ten  groin  compartments  of  the  eleven 
original  groins  were  not  filled  during  initial  construction,  natural 
longshore  transport  eventually  accomplished  this.  Over  the  years 
onshore-directed  wind-blown  sand  has  created  a  high  and  broad 
dune  that  is  fronted  by  a  wide  berm.  Beach  fill  accompanied 
construction  of  Groins  12-15.  In  December  1992,  a  northeaster 
opened  two  inlets  (called  Pikes  Inlet  and  Little  Pikes  Inlet)  directly 
west  of  the  groin  field.  Pikes  Inlet,  located  most  westward,  gradually 
shoaled  and  was  readily  closed  manually  (Bocamazo  and  Grosskopf, 
1999),  but  Little  Pikes  Inlet  enlarged  in  the  eroding  down-drift  area 
directly  west  of  Groin  15.  In  1993,  the  U.S.  Army  Engineer  District, 
New  York  (NYD)  closed  the  breach  by  hydraulic  fill  placement.  In 
1996-1997,  a  tapered  groin  transition  to  the  west  was  created  by 
shortening  Groins  14  and  15  and  adding  a  short  groin  (called  Groin 
14A)  between  them.  The  groin  tapering  together  with  beach  fill  has 
been  successful,  and  the  planned  renourishment  of  3  years  has  been 
extended  to  4  years. 

As  an  alternative  to  the  existing  condition,  the  NYD  considered 
shortening  of  groins  in  the  eastern  and  middle  portions  of  the  groin 


field  and  tapering  of  the  groins  on  the  western  end,  together  with 
continued  renourishment.  Groin  shortening  may  be  a  cost-effective 
solution  as  compared  to  beach  nourishment  from  an  external  source 
for  maintenance  of  the  beach  west  of  the  groin  field,  at  least  for  a 
certain  length  of  time  until  the  beach  and  dune  in  the  groin 
compartments  develop  equilibrium  with  the  modified  structures. 
Shortening  of  the  groins  has  the  potential  for  gradually  releasing  sand 
to  the  littoral  system  as  the  shoreline  and  dune  line  recede.  If  the 
beach  is  in  equilibrium  with  the  groin  field,  a  reduction  in  groin  length 
is  expected  to  lead  to  recession  of  the  shoreline  and  dune  as  the  beach 
is  reaching  a  new  equilibrium.  In  conclusion,  the  dune  and  berm  (as 
represented  by  shoreline  position)  are  expected  to  recede  if  the  groins 
are  shortened.  The  time  scale  of  this  process  is  not  known  a  priori  and 
is  an  output  from  the  application  of  the  new  model. 

Calibration.  The  model  was  set  to  match  the  annual  net  longshore 
transport  rate  at  the  site,  following  the  procedure  described  by  Kraus 
et  al.  (1994).  The  coefficients  for  sand  transport  to  and  from  the  dune 
were  set  based  on  the  experience  from  the  cross-shore  test  described 
in  the  previous  section.  Simultaneous  wave  and  water  level  data 
(from  the  Sand  Hook,  New  Jersey,  long-term  tide  station)  were  input 
to  the  model  to  calculate  long-time  dune  and  beach  berm  change. 
Simultaneous  waves  and  water  level  (including  storm  surges)  at  3- 
hour  interval  are  available  from  1  January  1980  to  31  December  1995. 
These  data  served  in  the  calibration  process,  for  which  the  measured 
shoreline  positions  of  1972  and  1983  were  used.  The  result  of  the 
calibration  is  shown  in  Fig.  10. 

Starting  with  the  measured  1972  shoreline  (short-dashed  line), 
the  main  features  of  the  measured  1983  shoreline  (long-dashed  line) 
were  reproduced  (solid  line).  In  particular,  the  continued  accretion  on 
the  up-drift  side  of  the  groin  field  is  replicated,  as  well  as  the  accretion 
in  the  eastern  and  central  sections  of  the  groin  field  and  the  erosion 
in  the  western  section.  The  model  also  predicted  the  amount  of  down- 
drift  erosion  well;  however,  it  is  located  too  close  to  the  groin  field.  The 
mean  calculated  transport  rate  alongshore  was  about  400,000  m3/yr, 
which  is  greater  than  the  long-term  average  estimated  value.  Overall, 
the  calibration  was  judged  to  be  satisfactory. 

Response  to  groin  shortening.  To  quantify  the  response  of  the  dune 
and  beach  to  groin  shortening,  the  model  was  run  for  20  years  for  two 
groin  length  alternatives,  one  with  the  original  lengths  as  described 
above  and  one  where  the  groin  lengths  were  reduced  by  30  m.  Both 
simulations  started  with  the  measured  1995  shoreline  and  included 
the  beach  fills  that  were  reported  to  be  added  in  1996-1997.  During 


Fig.  10.  Calibration  run  for  Westhampton  Beach.  Simulation  from  1972  to  1983. 
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Fig.  11.  Calculated  shoreline  change  at  Westhampton  Beach  over  20  years  with  existing  groins  and  with  reduced  groin  lengths. 


initial  construction  in  1996  and  1997,  about  2.7  million  m3  of  sand 
dredged  from  an  offshore  borrow  site  was  placed  from  Groin  7  to 
approximately  3000  m  down  drift  (west)  of  Groin  15.  For  the 
simulations,  the  fill  was  divided,  with  1.1  million  m3  within  the 
groin  field  and  1.6  million  m3  distributed  to  the  west  of  Groin  15.  The 
beach  section  consisted  of  a  berm  2.9  m  above  MSL  and  about  27  m 
wide. 

Fig.  11  (long-dashed  line)  shows  the  result  of  the  20-year 
simulation  with  the  original  groin  lengths.  As  expected,  accumulation 
continues  up  drift  of  the  groin  field  and  also  somewhat  inside  the 
western  part  of  the  groin  field.  Erosion  occurs  in  the  last  groin 
compartment  and  down  drift  of  the  groins.  The  result  of  the  second 
simulation,  with  the  shorter  groins,  is  shown  as  a  solid  line  in  Fig.  11. 
There  is  a  notable  difference  as  compared  to  the  results  for  the  original 
groin  lengths.  Shoreline  advance  (accretion)  up  drift  of  the  groin  field 
and  inside  the  western  side  of  the  groin  field  was  substantially 
reduced.  Likewise,  the  significant  recession  in  the  easternmost  two 
groin  compartments  and  down  drift  of  the  groin  field  was  also 
substantially  reduced.  The  reduced  down-drift  erosion  of  410,000  m3 
is  still  only  slightly  more  than  half  or  the  estimated  potential  of 
750,000  m3,  which  would  correspond  to  a  situation  where  the  entire 
beach  and  dune  system  would  move  back  the  same  distance  as  the 
groin  shortening,  i.e.,  30  m.  Thus,  it  is  estimated  that  the  sand  supply 
from  the  groin  length  reduction  will  persist  for  at  least  30  years. 

This  simulation  indicates  that  the  beach  and  dune  at  Westhampton 
Beach  will  recede  up  drift  and  inside  the  groin  field  and  that  the 
down-drift  beach  will  accrete  in  response  to  groin  shortening.  It  will, 
however,  take  several  decades  before  the  down-drift  beach  receives 
the  full  benefit  of  the  groin  length  reduction  and  before  a  new 
equilibrium  condition  is  established. 

5.  Conclusions 

A  new  numerical  model  has  been  presented  for  calculating  beach 
and  dune  change  under  combined  cross-shore  and  longshore 
transport  processes.  Predicting  this  change  involves  coupling  among 
the  dune,  berm,  and  shoreline,  and  it  describes  the  beach  as  a  balance 
of  dune  growth  by  wind-blown  sand,  dune  destruction  during  storms, 
change  in  berm  width  (which  enters  in  providing  the  surface  for  sand 
to  be  blown),  and  gradients  in  wave-induced  longshore  transport. 
Equilibrium  morphologic  concepts  were  included  to  improve  model 
predictive  reliability  and  stability  of  calculations  with  the  coupled 


non-linear  equations.  Simplification  of  the  governing  equations  and 
boundary  conditions  allowed  analytic  solutions  to  be  obtained  for 
examining  model  properties  and  performance. 

Sensitivity  tests  with  the  numerical  model  revealed  that  berm 
width  plays  a  central  and,  perhaps,  subtle  role  in  protecting  the  dune 
against  storms  through  its  storage  mechanism.  The  sensitivity  tests 
also  examined  the  role  of  increased  wave  height  in  storms  as  a 
possible  consequence  of  global  warming.  A  coastal  management 
conclusion  is  that  higher  volume  dunes  can  protect  the  coast  even  for 
higher  waves. 

Finally,  an  engineering  application  of  the  coupled  cross-shore  and 
longshore  transport  processes  model  was  made  for  evaluating  groin 
shortening.  Prior  to  development  of  this  model,  there  was  no 
technology  available  to  calculate  the  time  scale  of  combined  beach 
and  dune  change,  or  the  long-term  interaction  of  the  beach  and  dunes. 
A  conclusion  of  the  case  study  was  that  groin  shortening  would 
function  in  providing  the  desired  up-drift  feeder  beach,  and  the  sand 
would  be  released  over  several  decades  in  a  gradual  process  of 
adjusting  to  the  new  structure  lengths. 

This  study  has  brought  to  light  many  interesting  morphologic 
questions  and  focuses  them  within  a  quantitative  calculation 
framework.  For  example,  is  there  an  equilibrium  berm  width  for 
given  beach  sediment  grain  size,  transport  rate  by  wind,  and 
frequency  of  storms?  Such  questions  can  now  be  pursued  for 
improving  engineering  calculations  of  beach  and  dune  change  under 
realistic  conditions. 

Major  processes  presently  lacking  that  are  planned  to  be 
introduced  to  the  model  is  cross-shore  sub  aqueous  transport  by 
wave  action  and  by  wave  overtopping  and  overwash  of  the  dunes. 
These  further  developments  are  under  way  and  will  be  operational  in 
the  near  future. 

List  of  notation 

a  empirical  coefficient  in  runup  height  formula 

A  =1  Qo/Qwo 

B  =Dd/(Db  +  Dc) 

Cs  empirical  coefficient  in  dune  erosion  formula 

d50  median  grain  size 

Db  berm  height  (above  still-water  level) 

Dc  depth  of  closure 

Dd  dune  height  (above  berm) 

/  friction  coefficient  for  swash  flow 
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g  acceleration  due  to  gravity 

H0  deepwater  root-mean-square  wave  height 

kf  friction  coefficient  for  wave  traveling  over  the  berm 

I<w  empirical  coefficient  in  Bagnold  transport  equation  for  WBS 
L0  deepwater  wavelength 

q0  transport  from  dune  (erosion)  due  to  wave  impact  (offshore) 

qw  transport  by  wind  (onshore) 

qwo  maximum  transport  by  wind  (onshore) 

Q.  longshore  sediment  transport 

R  runup  height 

R'  adjusted  runup  height  with  respect  to  frictional  losses  over 

the  berm 

s  coordinate  originating  at  the  berm  crest  and  pointing 

onshore 

sb  =yB-yD 

t  time 

tf  qwot/w0DD 

T  swash  period  (taken  to  be  equal  to  wave  period) 

uc  front  speed  of  the  uprushing  wave  at  the  berm  crest 

u0  front  speed  of  wave  traveling  over  the  berm 

us  initial  front  speed  of  uprushing  wave  at  still-water  level 

u*  wind  shear  velocity 

VD  dune  volume 

w  =yB-yo  (berm  width) 

w'  =  w/w0 

we  equilibrium  berm  width 

w0  initial  berm  width 

y  shoreline  position 

yB  location  of  seaward  end  of  berm 

yD  location  of  dune  toe 

y50  distance  from  seaward  end  of  berm  to  where  the  WBS 

transport  has  reached  50%  of  its  maximum  value 
zD  dune  toe  elevation 

a  scale  coefficient  in  WBS  transport  formula 

a'  ow0 

A  change  in  a  quantity 

Ah  deviation  from  mean  sea  level  (includes  surge  and  tide) 

=  exp  (  —  aw) 

L  =exp(  — aw0) 


Subscripts 

o  related  to  wave  impact 

w  related  to  wind  transport 
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